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ABSTRACT 


A TDF rate integrating gyro in the strapdown mode is an attractive 
candidate for attitude sensing in a spacecraft navigation system, since 
it provides an additional axis of information for a relatively small 
increase in hardware complexity. A type of gyro which has not Been fully 
exploited is of the dry, tuned, TDF design, in which the spring constant 
of the suspension system is effectively cancelled by the "dynamic anti- 
spring" of a swiveling, rotating, gimbal. The use of this unconventional 
gyro in a digital rebalance loop is investigated. 
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CHAPTER I 


INTRODUCTION 

The use of gyroscopes as attitude sensors in spacecraft appli- 
cations is well known. Efforts to improve sensor performance while 
decreasing cost, weight, and power requirements are continuously 

being made. One such effort which has been successfully exploited 

1 

is the strapdown concept, in which the gimballed platform upon 
which the gyros reside is functionally replaced by an electronic 
rebalance 'loop for each gyro. This loop senses float motion away 
from a null position with respect to the case and applies the correct 
amount' of current to the torquer coils to drive the float back t<3 
null, or rebalance it. Since the gyro floats’ attitude with. respect 
to the case is invariant, the case may be mounted to the- spacecraft 
frame, eliminating the need of a separate platform. Thus the gim- 
balled platform is said to be replaced by a computer, or analytical 
platform, and the quantity of sensor, and spacecraft, motion is found 
by measuring the torquing current. The advantage of the strapdown < - 
concept is that mechanical components are replaced with. electronic 
components with an Improvement in cost, weight, and reliability. 

Although single-degree of freedom (SDF) floated gyros are pre- 
sently used' in strapdown systems, two-degree of freedom design would 
give another axis of information at a nominal increase in cost and 
weight. Conventional TDF floated gyros for platform use have not been 
adapted to- strapdown applications because of difficulty in torquer 
design. The unconventional mechanical design described in this report 
utilizes a tuned suspension which eliminates the need for rotor 
flotation, high quality bearings, and elaborate gimbals. This design 
utilizes an inside-out construction in which the gyro rotor is ex- 
ternal to its support and drive mechanisms, making torquer design easy 
to implement. In fact, this dry, tuned TDF gyro is ideally suited to 
strapdown applications since its limited angular range requires a 
rebalance loop. 
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In rebalance loops for strapdovn applications, a digital tor- 
quing scheme commonly used is known as a pulse torque servo amplifier 
(PTSA). Accurately known current pulses are delivered to the tor- 
quer coil. The advantages of using a PTSA are twofold: 

1. Because the current is either zero, or a positive 
or negative maximum, non-linearity of the torquer 
coil is eliminated as a source of scale factor 
error. 

2. ' The current pulses may be counted to indicate the 

restoring torque necessary for rebalance, and con- 
sequently the sensor motion. This technique elim- 
inates the need for additional analog to digital 
conversion, provided the pulse quantization is fine 
enough for the required data resolution. 

Toward the goal of using a TDF dry tuned gyro as a strapdown 
instrument, a feasibility study for a PTSA loop is included in this 
report. For a better understanding of the control problem the trans- 
fer function of a tuned gyro is derived. As an example of a practical 
rebalance loop, the U„T, binary width modulated loop can be modified for 
use on a TDF dry tuned gyro. 
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DYNAMICS OF A DRY TUNED TDF GYRO 

I. DESCRIPTION OF THE TUNED GYRO 

In the early 1960's the need for a relatively inexpensive, 

light' weight, rate-integrating gyro was recognized for use ‘as an in- 
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ertial sensor. 5 5 5 Conventional floated gyro technology had be- 
come sophisticated to the extent that additional improvements were 
very expensive. Undesirable effects, such as rotor mass unbalance 
caused by asymmetry of the windings on the spin motor, spring restoring 
torques caused by the motor power leads, and temperature' effects caused 
by thermal gradients of the floatation fluid, could not be reduced 
easily and cheaply. Thus the stage was set for a radical departure 
from conventional gyro design philosophy. 

This new philosophy produced a "tuned" gyro (Figure 2-1A). in 
which the rotor, coupled by a rotating, bearingless, suspension system 
to the drive shaft and motor, is "outside" its support mechanism, 
rather than inside, as in conventional design. The suspension system 
consists of a single ring, or gimbal, connected to the shaft and rotor 
by elastic restraints which have finite stiffness in torsion (twist) 
and Infinite stiffness in flexure (bending). The axes of these re- 
straints are orthogonal. It is the inertial reaction torque generated 
by the motion of the gimbal relative to the rotor and shaft effectively 
cancelling the elastic torque generated by twisting the restraints, 
which allows the rotor to be free, that is, torsionally decoupled from 
the shaft about its axes of freedom. This inertial reaction torque is 
sometimes called a dynamic antispring. 

The tuned gyro is termed unconventional, or inside out, to dis- 
tinguish it from the conventional gyro (Figure 2-1B), in which the 
rotor, and its spin motor, to which it is directly coupled, are inside 
the support mechanism. These gimbals do not spin with the rotor, as 
in the tuned case, rather they decouple the rotor from the case through 
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high quality hearings mounted on orthogonal axes, A gyro used in 
navigation and guidance systems must have a highly refined suspension 
for the rotor. One widely used means of "floating" the rotor con- 
sists of suspending a sealed housing containing the rotor and a spin 
motor in a high density fluid. The tuned gyro, which does not require 
flotation, is called "dry". 

The principle difference between conventional and tuned gyros 
is the physical phenomena utilized. Conventional gyro design has 
exploited every means to provide the. perfect hearing, while tuned 
gyros use the dynamics of the suspension system to give decoupling. 

II. THE DYNAMIC ANTISPRING 

The glmbal motion which generates the inertial reaction torque 
is simple harmonic with a frequency equal to twice the spin frequency. 
Figure 2-2 illustrates this motion, in which the gimhal plane ex- 
periences a complete period of its oscillation during one-half a rotor 
revolution because the axes of the restraints coincide after one-half 
revolution. 

Analytically, the antispring may he observed in the equations 
of motion resulting from the gimhal action described in the preceeding 
paragraph. These equations of motion .are • ' 

(A + A /2) <L + DO + [K - N 2 U - C jm 0 Y + 
g X X g S a 

a(c ♦ A g > + T D e Y - -|[A g e x + a g ue Y - 

(2A - C ) N 2 0 v ] cos 2N - |[A 0 - 2AN© - 

ggY 2gx g A 

(2A - C ) N 2 0 v ] sin 2N + G v 

g g Y x 


and 


^(2-D 
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Figure 2-2. ' Simple Harmonic Motion of the Gimbal 
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(A + A g /2) 0 Y + D0 y + [K - N 2 (A g - C / 2)] 0 y - 

B(C + A g ) « x - T d S x - |[A g S r - 2A g H0 x - 

(2A - C ) W 2 0„] cos 2N - i[A 0__ + 2A N0„ - 

g g Y J 2 g X gY 

(2A - C ) N 2 0„] sin 2N + G„. 
g g A X 


I (2-1) 

| continued 


© Y and Qy are angular deflections of the rotor about the case-fixed X 

and Y axes, in response to G and G , externally applied torques to 

rotor about these same axes. H is the spin frequency of the gyro, and 

all other terms are physical constants; C, A, C g , and A g spin and cross 

axis moments of inertia of the rotor and gimbal, respectively; D, 

viscous damping; K, spring constant; and T , the rotor-to-case drag 

w 2 

coefficient. The dynamic antispring .term, - N (A - C/2), is always 

&> S 

negative because A must be greater than C /2 (except for a gimbal of 
& & 

zero thickness, where A = C /2). If the spring constant equals the 

& s 

dynamic anti spring, the torque- free operation is obtained provided the 
second harmonic terms in (2-l) can be neglected, and the gyro is said 
to be tuned. If these terms cannot be neglected, as in the case where 
Gjj. or Gy contains 2N frequency components, then an unstable steady- 
state solution is obtained, because there is no positive spring con- 
stant term to cancel the antispring term in these second harmonic 
coefficients. This phenomenon is called the 21 sensitivity of the 
tuned gyro. 


III. TRANSFER FUNCTION OF A SINGLE GIMBAL TUNED GYRO 

A useful mode of operation of the gyro as a sensor is the strap- 
down mode, in which the gyro is strapped directly to the vehicle, 
rather than being mounted on a gimballed platform, and a rebalance 
torque applied to the rotor to maintain its attitude with respect to 
the vehicle fram. Attitude sensing is obtained by measuring the torque 
current necessary to rebalance the gyro. Accurate attitude information 
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depends on knowledge of the gyro transfer function. To give the 
reader a better understanding of the dynamics involved in a tuned 
gyro, a detailed derivation of its transfer function will be 
given . 

In this section the coupled differential equations of motion 
in rotating coordinates will be derived for a single gimbal gyro. 

These differential .equations will then be transformed into a trans- 
fer function relating case-fixed variables via a method using sym- 

6 T 

metrical components and a complex coordinate system. 5 The moment 
equations are derived in a shaft-fixed frame rotating at spin velocity 
because the dynamics of the rotating suspension system are easier to 
visualize in this frame. The method of symmetrical components provides 
a useful means to analyze systems with asymmetries, while the complex 
coordinate system is easy to use for the transformation between a 
rotating and stationary frame. 

Moment Equations 

Four coordinate systems (Figure 2-3) are needed for this de- 
rivation: 

(X,Y,Z) A case-fixed system, with its Z axis along 

the shaft spin-axis. 

(x,y,z) A shaft-fixed system, with z along the shaft 

spin-axis and x along the inner restraint. 

(x ,y ,z ) - A gimbal-fixed system, with x along the inner 

® ® ® restraint axis and y the outer. 

g 

(x' ,y',z') - A rotor fixed system, with y" along the outer 
restraint axis and z" the rotor spin axis. 

Figure 2-4 shows the detailed relation between case and shaft 
systems, in which the x and y axes rotate, at N, the spin velocity, 
in the plane formed by the X and T axes. The relation between the 
absolute angular velocity of the shaft, resolved along the shaft 
set, and the absolute case velocity, resolved along the case set, 
is expressed with the aid of the Euler-angle transformation as 



X, Y, Z - Imbedded in Case (non-rotating) 

x, y, z - Imbedded in Shaft (rotating at K radians /sec) 

x , y z - Imbedded in Gimbal (rotates and moves with gimbal) 

£ &> £ 

x", y", z' - Imbedded in Rotor (rotates and moves with rotor) 


Figure 2-3. Coordinate Frames Defined 



Gimbal 

Shaft 



Figure 2~k. Relating Case and Shaft Coordinate Frames 
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In Figure 2-5, the rotor is slightly tilted with’ respect to the 
shaft, such that angular displacements 0 and 6 "* have been made about 

•fl- y 

the x and y' axes, respectively. The absolute rotor angular ' velocity , 
of*, resolved along the rotor set, may be expressed as 


absolute rotor 
velocity along 
rotor set 


absolute shaft 
velocity along 
rotor set 


rotor velocity 
+ with respect 

to shaft along 
rotor set 


or 


w 


-(r ) 


(0 


rs(r) * 


(2-3) 


The term ma y expressed in terms of as 


— (r) 


cos Q ' 0 sin 0 

y y 


sin 0 " 0 cos 0 * 

y y 


0 • cos 0 sin 0 
x x 

0 -sin 0 cos 0 
x x 


w 


x 


tn. 




(2-4) 


Xn a tuned gyro, the angular deflections are small (in practice less 
than 5°) , since this is a null type instrument, in which rebalance 
torques are used to make 0^ and 0^ approach zero. Therefore, the 
following small angle approximations are valid: 
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Figure 2-5. Relating the Rotor and Shaft Frames 
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COS 6 ~ 1, COS 0 - ~ 1 

X y 

sin 8 ~ 8 , sin 0 " ~ sin 8 
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(2-6) 


The second term on the right side of (2-3) allows time variations of 

8 and 8 ' and is 
x y 


— — 


[- 


— - 

cos 0 ' 


cos 90 


0 

y 




X 

cos 90 

0 + 

cos 8 

0 
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X 

X 

y 

y 

» • 

sin 8 ' 


sin 0 


0 0 +00 

y 


X 


y x x y 


(2-7) 


using the small angle approximations. Adding (2-7) to (2-6), and 
observing that 


0) 

z 


N » 


0 9 
y x 


0 .e 
X y 


w 0 

y x 
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(2-3) becomes 



r— - 


— 

4 — 


"V 


to — N8 + 
x y 

0 

X 

to" = 

to- 

y 


to + W0 + 

y x 

8 

y 


& 

(0 .. 
z 


N 



(2-8) 


Figure 2-6 shows a top view of the gyro and a detailed relation be- 
tween the shaft, gimbal, and rotor sets. Note that if the rotor is 

not tilted, 8 and 0 - equal zero, and these sets are coincident. 

5 x y ^ 5 
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The absolute angular rate of the gimbal , oj resolved along the 
gimbal set, is 


Absolute Gimbal Absolute Shaft Gimbal Velocity 

Velocity Along = Velocity Along + "with respect to 
Gimbal Set Shaft Set shaft along 

Gimbal Set 

or 


w 

-g 


^g) 


+ 


^gs(g)' 


(2-9) 


Using the same rationale-as previously. 



i o -e 


0) 

^g) - 

yg 

o i e 

xg 


X 
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e -e i 


(0 


yg xg 


z 


and 


^gs(g) 


0 

xg 

0 


e 

zg 


( 2 - 10 ) 


( 2 - 11 ) 


The inner elastic restraints cause the shaft and gimbal axes, x and 

x , to coincide, while the outer restraints cause the gimbal and 

rotor axes-, y and y", to coincide. Thus the angles 

g 


8 a 0 , 

Xg X 


6 

yg 


- o, 


and their derivatives. 
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e = e , 

xg x’ 

e =o 
yg 


may be used -with (2-10) and (2-11) to give 


co 

“g 


— — 
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CO'*'* 


CO +0 

Xg 


X X 

hi" 
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CD 

& 

1 

3 

yg 


y • x 

CO'" 
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Zg 




( 2 - 12 ) 


Figure 2-JA shows a free body diagram of the gimbal used to 
derive its moment equations with respect to the gimbal coordinates. 
These moment equations, about the x and y axes, respectively, are 


and 


T - T 
xg xg 


' - A !£," 
g 


+ (C - B )■ co" co' 

xg g g yg 


zg 


T - T " = B co' 

yg yg g. 


' - (C - A ) hi" to' 

yg g g xg 


Zg 


>(2-13) 


where the T and T are from the rotor via the outer restraints, and 
xy yg 

T > and T ^ from the shaft via the inner restraints , and A , B , and 
x g y g g g 

C represent moments of inertia about the x , y , and z axes, respec-.' 
g g g g 

tively. A and B are called cross-axis, while C is called a spin- 
8 g g 

axis moment of inertia. The twisting moments yield the following 
torque equations. 


T . = K 0 „ + D 8 „ 

x g xx xx 


T =K6 , + D6 , 

yg ,y y y y 


!> (2-lU) 


where 0,, = 0,0., s: 0 > and K and K are spring constants, and D 
x x y y x y x 

and D viscous damping. Substituting (2-lU) into (2-13) gives 

v 
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T -Am" + (C - B ) m" m" + K 9 + D 0 

xg y xg g g yg zg xx x x 


\> (2-15)' 


-T „ = B &-* - (C - A ) u" a" - K 0 - D 9 

yg yyg g g xy zg yj yy 


From Figure 2-JB, the moment equations of the rotor about the- 
x' and y" axes are 


M , — Aco' 1 ,. + (C — B) in' ^ m' ^ + T^ + T_ sin 9 
x x y • z . x • D y 


and 


M , = Bai' . - (C - A) u' . bi' „ + T „ - TL sin 9 , 
y y x z y . D x 5 


(2-16) 


where and are externally applied moments, A, B, and C, rotor 

moments of inertia, ! „ and T ^ torques applied through the sus- 

x y 

pension, and the rotor drag torque. For small 0^ and , 


M - M and M ,, - M , 
x - x y y 5 


T ^ - T and T „ - T , 
x xg . y yg 


and (2-l6) becomes 


1 ' 

M x = Aa) V + (c ~ B) “V“V * V " x s + (c g " V z § + 


K6 + D 0 + T_0 

xx xx D y 


M = Bm' „ - (C - A) oi' .to' „ + K 9 + D 6 - T_8 , 

y y x z yyyyDx 


j>(2- 


17) 


where the first of (2-15) has been substituted for T , the second of 

xg 

( 2—1 4 ) for T , and the small angle approximation used for sin 0^ and 

yg. 

sin 0 ... Substituting (2-8) and (2-12) into (2-17) to get all angular 

y 

velocities in terms of shaft rates gives 
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M = A(a) + 0 - N0 ) + (C - B) (os +9 + W6 ) N + 

x xx y y y x 

A (os + 0 ) + (C — B ) (co + N8 ) I + K 0 
gx x ' g g y x ' xx 


+ D 0 + T,0 , 

xx D y 


M = B(o) +0 + N0 ) - (C - A) (u H- 0 - N0 ) N + 

y y y x' x x y 

K 0 + D 0 -10 . 

y y y y D x 


? (2-18) 


collecting terms with angles, -velocities, and their derivatives as 
common factors, gives the coupled differential equations of motion of 
the rotor in shaft fixed coordinates: 


0 .(A + A ) + D [K + It (C - B+ C - B )] + 

X g xx ^x x g g 

• ♦ 

0 K(C - B - A) + T +0 (A+A ) + 0N(C-B+ C - B ) 
y yD w x' g v y v g y 


= M 


x 


>(2-19) 


0 (B) H D + 0 [K + N (C - A)] + 

y y y y y 

• • - 

0 N(C — B - A) - 0 T + w B - o> N(C - A) 
x x D y x 


= M 


y 


Transformation Using the Complex Method 

To transform the coupled differential equations of (2-19) to a 
transfer function relating variables in a case-fixed frame, a complex 
plane technique will be used. First, define the following symmetrical 
components: 

I *» (A + B + A )/2 

g 

I' = (A - B + A )/2 
g 
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D * (D x + D y )/2 

D" «* (D - D )/2 
x y" 

K = (K + K )/2 
x y 

K" = (K - K )/2 
x y 

J * (2C - B - A + C - B )/2 

g g 

J' = (A - B + C - B )/2 

g g 

L = C - B - A 
L" = 0 


/ (2-20) 


P = (A + B + A )/2 
g 

P' = (A - B + A )/2 
g 


R = (2C - B - A + C - B )/2 

g g 

R" = (A - B + C - B )/2 
g g 


Substituting (2-20) into (2-19) 


M 


M 


y 


0 (l + I") + 6 (D + D") + 8 [K + K" + N 2 (J + J**)] 
x x x 

• • 

8 N(L + L")0 XL + 8 (P + P") + 8 H(R + B') 
y y D x y 

0 (l - I") + 0 (D - D-) +0 [K - K- + U 2 (j - J")] 

y y y 

• * 

0 N(L - L") - 0 T + w (P - P") - w N(R - R") . 
x x D y x 


>( 2 - 21 ) 


Next the equations of (2-2l) are combined into a single equation 

• 1 1 > 

by multiplying the second by j - / -1 and adding it to the first, 
using the following complex relations: 
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their conjugates. 
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xy x 
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and derivatives, and conjugate 'derivatives , 
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xy 
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M - 10 + I"e + (D - jKL) 6 + (D" + jML') 9 

xy xy xy xy xy 

(K - h 2 J - jT ) 0 + (K" - W 2 J') 8 + Pto + 

° D xy xy xy 
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xy xy xy 


/ (2-23) 


The. Euler angle transformation from a shaft-fixed (rotating) 
frame to a case-fixed (non-rotating) frame for- angles, rates, and 
their derivatives, are 



0 

. xy 
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xy 


= 0 XI E 


-JHt 


‘ <®xr - J”y ^ 



(2-2k) 
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- 1 -jift 

y f 

«_ - (+XY - OTv>' ,,Nt ' 


xy 


^XY 


while conjugation of the above relations simply conjugates each term of 
each relation, which changes the sign of the imaginary component. The 
externally applied moment minus the rotor damping is 

i l 

V " - VxF^- < 2 - 2 5> 

Using the relations of (2-2U) and (2-25) in (2-23) gives 

{itG^ - j2H0 xy - W 2 ©^] + (D - jKL) (Gjjy - jwe^) 

+ (K + N 2 J - jT D ) + 

{rtl x!c + - At* 3 + (D ' + 3HL ' ) (t xY + jb ®xy) 

+ (K" + N 2 j") &_} e +JWt = ' (2-26) 

AX 

- P(i xy - P^) - jNE^xy) e" JWt + 

{-p' (tjy + jw?^) - e +JNt 

which is the complex form of the rotor moment equation in case-fixed 
coordinates. The following Laplace transformations are useful: 


L[f(t)] * F(s) 

L[f(t)e +j21ft ] = F(s - J2N) 

L[f (t )] - sF(s) l 

L[f(t)3 = s 2 F(s) 

L[f(t)e J2Ht ] = (s - J2H) F(s - j23?) 

L[?(t)e j2Mt ] = (s 2 - jPs - P 2 ) F-(s - j2N) 

where initial conditions are assumed zero. Equation (2-26) is multi- 

iNt 

plied through by e d and Laplace transformed using (2-27), giving 
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U(s 2 - J2Ns - W 2 ) + (D - JHL) (s - JH) - JT R + K + N 2 J + 
D r s } ©^(s) + {I'(s 2 - J2Ns - N 2 ) + (D" + JHE) (s - JH) 

K" + hV} 0^,(3 - J2H) = M^U) - {P(s - JH) - JHR} <^ (s) 
- {P'(s - JH) + JHR'} (s - J2H). 


The coefficients of ©^s, ©^(s - j2N) , - ^(s), and - ^(s - J2N) 
are defined as 


Z (s) = I(s 2 - J2Hs - H 2 ) + (D - JHL) (s - JH) + K + N 2 J - JT R + D R s 
Z (s) = I'Cs 2 - J2Ns - H 2 ) + (D' + JNL") (s - JH) + K' + U 2 J 
Z 3 (s) = + {P(s - JH) - JMR} (2-29) 

Z 4 (s) = {P"(s - JH) + JHR"} 

Substituting (2-29) into (2-28) yields 


Z 1 (s)0 XY (s) + Z 2 (s) ©^(s - J2N) * M^s) - Z 3 (s) ^(s) 

"V s) ^XY (s “ j2K) * (2 " 30) 


Equation (2-30) is a statement in the complex frequency domain relating 

complex rotor angle 0 (s), and its conjugate at twice spin frequency, 

_ XX « 

©•^(s - J2H), to the complex case input rate, ^^(s), and its 2H con- 
jugate 4>jry( s - J2W), and externally applied complex moments, M^fs). 
This equation may be solved for either the direct or 2H frequency com- 

g 

ponent by the following means 0 : 

1. Conjugate (2-30), using a dummy variable for s. 


2. Substitute s - J2H for the dummy variable, 

3. The new equation, which is 


z i (s - J 2H) V s - ^ + V s - J 2 ”> ®xr = ( 2 - 31 ) 

* « 

^(s - J2N) - Z 3 (s - J2H) ^(s - J2H) - Z^s - J2H) ^(s.) 



2k 


may "be solved simultaneously with (2-30 ) to give either 
Qjry(s) or 0^(s - j2N). This solution , in matrix 


form 

S XY (s) 

, is 

Z 1 (s) Z 2 (s) 

-1 

“HC(a) 

9^(3 - J2H) 


Z 2 (s - j2H) Z^(s - J2N) 

4» 

\y(s - J2N) 

— 


— - 


— _ 


(2-32) 


V s ) ■ ■ 

Z 2 (s) 

-1 

« 

Z 3 (S) 

Z 4 (s) 

0 

*xr (s) 

• 

Z 2 (s - J2N) 

Z x (s - J2B) 


Z 4 (s - 32H) 

Z 3 (s - o2K) 


^y(s - <J2N) 


Justification for the preceeding operation is given in Appendix A. 

Solution of (2-32) is simplified using the following approx- 
imations : 


1. Rotor to case drag, and damping coefficients are ne- 
gligably small, therefore 

T D = D = = D R = 0. 

2. By symmetry, 

A a B, 


and K = K 

x y 

therefore, (2-20) becomes 

I = A + A /2 
g 

= A g /2 
K' - 0 
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J = C - A + 1/2 C - 1/2 A 

g g 

J" = C J2 - A /2 
g g 


L = C - 2A 
L" = 0 


P = 2A + A 
■g 



R=C-A+C/2-A/2 

* g g 

R" = C /2 - A /2 
g g 


and ( 2 - 29 ) may "be simplified to 


> (2-33) 


2. (s) = s 2 (A + A /2) - jWs(C + A ) + W 2 (A - C /2) 

1 g g g g 

Z (s) = s 2 A /2 - jNsA - N 2 (A - 0 / 2 ) 

2 g g g g 

Z (s) = s(A + A /2) - jrr(C + C /2) 

^ s § 

Z, (s) = sA /2 - jlf(A - 0 / 2 ). 

4 g g g 

conjugation and substitution of s - j2N in (2-3*0 yields 


>(2-3*0 


Z-, (s - J2R) = s 2 (A + A /2) - jNs(4A - C + A ) + 

1 g g 

K - N 2 (i*A - 2C + A -0/2) 

g g 

Z 2 (s - j2H) = s 2 A g /2 - jRsAg - N 2 (A g - C g /2) 

Z_(s - j2K) =s(A+A/2) - j I (2A-C + A - C /2) 
b g g g 

Z, (s - J2N) = sA /2 - m /2. 

' o o 


>(2-35) 


Solution of (2-32) for ©^(s) gives 

®XY (s) = Z l (s “ j2Il) M XY (s) " Z 2 (s) \y (s “ j2N) ^ 

- [Z x (s - J2H) Z 3 (s) - Z 2 (b) Z 4 (s - j2N)} ^(s) 

- [Z r (s - j2H) Z 4 (s) - Z 2 (s) Z <s - J2N)] ^(s - J2N) 

[Z 1 (s) Z ± (s - j2N) - Z 2 (s) Z 2 (s - j 2M ) ] , J> (2-36) 
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which is the complex form of the output angle, in case-fixed coordin- 
ates, as a function of input rates and moments. 

Transfer Function 


Reduction of (2-36) to a useful transfer function form, requires 
the following assumptions: 

1. The rotor is thin, that is,.C_= 2A, which allows 
simplification of the terms Z.. (s - j2R) and 
Z 3 (s - J2N). • 

2. ' Gdmbal inertias are negligably small in comparison 

with rotor inertias, which allows the following 
simplifications : 

A. The product Z^Cs) • Zg(s - j2N) in the 
denominator may be neglected. 

B. The effect of M^.(s - j2N) may be neglected 
because Z^(s) is small. 

C. The product Zg(s) • Z^(s - j2W) in the co- 
efficient of ^^(s) may be neglected. 

D. Because of A. , the coefficient of ( s - j2N) 
may be reduced to 

V s5 ) V s ) 

“ V's7 + Z 1 (s) • (s - j2N) 5 


which, since it contains only gimbal terms 
in the numerator, may be neglected. 

Use of the preceeding assumptions and the condition of tuning, 

K = U 2 (A - C /2), (2-37) 

S S 


to eliminate the constant terms in Z^(s) and Z^(s j2Sf), reduces (2-36) 
t'o 

, . Uy^(s) 

0 XY U ' As(s - j2R) " 


s 


(2-38) 
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Equation (2-38) may be transformed using (2-22) into 


0 x (s) 


M^sVA 

~~p Q 

s + hu 


(2N/A) • My(s) 
s(s 2 + In 2 ) 


¥ s) 


®y(s) 


M^(s) * 2K/A 
s(s 2 + 4 h 2 ) 


^(s)/A 
+ — — 
, s + 1*N 


4> y (s) 

s 


> (2-39) 


7 

which is the transfer function of a two-degree-of-freedom rate in- 
tegrating gyro. . 

Transfer Function Including 2N Rate and Moment Inputs 


Some of the preceeding approximations used to simplify the 
transfer function of a tuned gyro are inappropriate. In a gyro used . 
in the strapdown mode. Figure 2-8, the rebalance torque is applied 
with a magnetic force, and proper design must provide an efficient 
magnetic circuit for the torquer flux. Consequently,' the rotor is 
not thin, that is, 2A =/ C. For the rotor of Figure 2-8, 2A is about 
10 % greater than C, and the quantity 2A - C cannot be neglected in 

Z.(s - j2N) and Z (s - j2W), since it is large compared with A - C /2. 

-L . j e> S 

Although higher ordered terms of Z 2 (s) • Z^Cs - j 21T ) may be 
neglected in comparison with those of Z^(s) • Z^(s - j2N), its constant 
term must be considered. The condition of .tuning, which forces the 
constant term of the characteristic equation (Z^(s) Z^(s - ,5211) - 
Z 2 (s) • Zg(s - j2N)) to zero is 


K n U 2 { (A 

g 


C /2 + 2A.-« 
g 


C) -t [(A - C /2) 2 
g g 


P 1/2 

+ (2A - c) 2 ]- ). 

(2-%)) 


Other terms of the characteristic equation are 


A 2 s^ - jJT4A 2 s 3 


H 2 hA 2 s 2 + «)Un 3 A 


[(A g"V 2) + ( ’ 2A 


c) 2 ] 


1/2 


(2-4l) 
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in -which the constant term is negligably small, leaving 

p p p 

A s (s - j2N) 


(2-42) 


as the characteristic equation. The numerator of the transfer func- 
tion is 

-A 2 s[s 2 - j4us 4 n 2 ] * (s) + As (s - J2N) M^s) 

+ JNAU - C 72) s(s - jN) L(s' - J2H) (2-43) 

g g XX 

+(A /2) s(s - j2N) M^s - J2H), 

■where the constant terms are considered negligably small, because they 

contain inertial differences due to finite thickness, i.e., A - C /2 

•_ g g 

and A - C/2, or in the case of the ({^(s - j2N) term, the product of 
these differences. The transfer function may be resolved into its X 
and Y components, giving 


V. s > - - 


<j> x (s) sM^/s) - 2N • My(s) 

s As(s 2 + 4U 2 ) 


A / 2 
-Jg.— 


N(A - C /2) 
g g 


sM^s - J2N) + 2W • My(s - J2N) 

As(s 2 + 4N 2 ) 

(3N s 2 + 4N 3 ) <^ x (s - J2H) - s 3 <j> y (s - J2N) 
As(s 2 + 4K 2 ) 2 


and 


>(2-44) 
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<iy(s) 

Vs) + 


2N 


* M^s) + sMy(s)- 
As(s 2 + 4 n 2 ) 


A /2 21 • M^(s - j2K) - s • My(s - j2N) 

A As(s 2 + 4N 2 ) 

N(A - C J2) ~s 3 <f> (s-- J2K) - (3Ns 2 + 1+1I 3 ) <f> (s - j2N) 

e s . A y. 

As(s 2 + 4 n 2 ) 2 


Equations .(2-44) contain all the terms of (2-39) plus additional terms 
to account for external moments and input rates at twice' the spin freq- 
quncy. Of particular importance are the terms related to 2N rate in- 
puts. An impulse rate input about either axis at 2N (a sinusoidal 
angular displacement of constant amplitude and frequency of 2N) gives 
outputs about both axes which have a non-zero average value. This 
phenomenon is known as 2N rectification, and generates a drift error 
in the gyro used as a sensor. The terms related to 2N moment inputs 
are not particularly important, because moment inputs are controlled 
by currents in the torques, and the designer can eliminate currents 
at this frequency. A block diagram of (2-40) is shown in Figure 2-9. 

IV. TRANSFER FUNCTION OF A TWO GIMBAL TUNED GYRO 

For the single gimbal gyro with a thick rotor (2A > C), the 
tuning condition (2-40) contains rotor inertias, indicating that rotor 
dynamics enters the design of a single gimbal tuned gyro. The addition 
of a second gimbal and pair of restraints (Figure 2-10) allows the 
designer to eliminate rotor terms from the tuning equation and 2N 
rectification from the output. 

Moment Equations 

The second gimbal adds another coordinate system to the four 
considered for the single gimbal case. 




Figure 2-9. Block Diagram, of .Tuned Gyro Dynamics' - Including 2N Components h 
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(x , y , z ) Fixed in the outer gimbal, with 

s s s x p along the outer restraint 

axis and y „ along the gimbal 
to gimbal ® restraint , 

also, the system for the inner gimbal is now designated by 





Fixed in the inner gimbal, with 
x - along the inner restraint and 
a l° n S the gimbal to gimbal 
r§straint . 


Equations (2-2) and (2-8) for shaft and rotor .velocities are valid 
for the two gimbal gyro . Following the development of ( 2-12 ) , for 
the inner gimbal. 




Referring to Figure 2-11, if the restraints from gimbal two to the 

rotor have the same restoring torque as those between gimbal one and 

the shaft, the angular deflection of the gimbals about the x-axis is 

! 

one-half that of the rotor. About the y-axis, the inner gimbal has 
no deflection, while the outer experiences full rotor deflection. 
Consequently, for the inner gimbal: 



3b 



C . Front 


Figure 2-11. Relative Deflections of Gimbals and Rotor 
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.xgl 

5 Xgl 


0 /2 
x 


8 /2 
X 


e 


= e 


ygi ygi 


= o 


0 


zgl 


~ o 


> (2-47) 


and (2-45) becomes 


~gl 




“ + ® /2 

xgl 


X X 



a + K0 /2 

ygi 


y x 

“zgl 


N 


(2-48) 


while for the outer gimbal, 


xg2 


xg2 

3 yg2 


0 / 2 
x 


0 / 2 
x 


= 0 


yg2 


= 0 


zg2 


y 

.o , 


> (2-49) 


and (2-46) becomes 


% 




w - N0 +0/2 

xg2 


x y x 

u) ^ 

— 

to + N0 /2 + 0 

•yg2 


y x y 

CO ^ 


u 

zg2 

w f i 




(2-50) 


Figure 2-12 A, B, and C shows free-body-diagrams of the inner 
gimbal, outer gimbal, and rotor, respectively, of the two gimbal gyro,- 
used to derive the moment equations of the rotor. The moments t>f 
inertia about the respective x, y, and z axes are designated A, B, and 
C, while the subscripts 1 and 2 denote inner and outer gimbals, with 
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no subscript for the rotor. The inner and outermost restraint torsional 
spring constants are and K , vhile that of the middle restraint is 
K . Viscous damping coefficients are defined in a similar manner. 

From Figure 2-12a, the moment equation about the x ^ axis is 


M 


xl 


- A l“xgl + ^ C 1 B l^ “ygl^zgl + K xl 0 xgl' + I) xl 9 xgl ’ 


ygl zgl 


xgl 


-(2-51) ' 


where is the moment delivered from outer to inner gimbal via the 
middle restraints, which have zero bending. Consequently,- for the outer 
gimbal. Figure 2-12b, the moment equations about x 0 and y 0 are 


and 


” M xl A 2 t0 xg2 + ^ C 2 “ B 2^ t0 yg2 t0 zg2 + K x2' 9 xg2 ” 0 x ^ 

o • 

+ 'D _(0 0 - “0 ") , 

x2 xg2 x 


> ( 2 - 52 ) 


M y2 B 2 a) yg2 ~ ^ C 2 A 2^“xg2“zg2 + K y^ 9 yg2 " 0 ygl^ 

+ D y^ 9 yg2 ~ 9 ygl^ ’ 


respectively, where is a bending moment from the rotor through 
the outer restraints. For the rotor, from Figure 2~12c,- the moment 
equations are 


and 


M ** - Aw" „ + (C - B)w" ,tu" „ + K _(e * - 0. J 
x x • y z x2 x xg2 

ft * 

+D O (0 „ - 0 0 )+ T 'sin 0 , 

x2' x xg2 . D y 


>(2-53) 


M , - M _ = BuT ^ - (C - A) to" ,<*>' , 

y y2 y x z 


bout the .x - - and y* axes. Equations (2-5l) and (2-52) may be combined 
ith (2—53) to give 
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M „ = Ago'' . + (C - B)to' .to" . + A to + (C - B )to O to _ 

x x y z 2 xg2 2 2 yg2 zg2 

• « 

+ A 1 oj + (C - B )to go + K 0 , t D 9 

VL xgl 1 1 ygl zgl xl xgl xl xgl l (2-5U) 

M . = Bio' . - (C - Ajco" .to" „ + B-io - (C 0 - A )io O to 0 

y y x z . 2 yg2 2 2' xg2 zg2 

+ K(6 . - 0 ) + D ( 0 -9 _ ) . 

y yg2 ygl y yg2 ygl 

Using ( 2-47 ) and (2—49 ) , and the small angle approximation which allows 

M . ~ M and M , : Mi (2-52) "becomes 
xx y y 


M = Ato' „ + (C - B)io' .to' . + A to _ + (c. - B n )<o _u _ 
x x y z 1 xgl 1 1 ygl zgl 

♦ • 

+ A 0 co _ + (C 0 - B )to „(o 0 + K .. 0 /2 + D-.0 /2 

2 xg2 2 2 yg2 zg2 xl x xl x ( 2 - 55 ) 

• • 

M = Bio' . - (C - A)to' ,10' . + B„to * - (CL - A^)io 10 _ 
y y x z 2 yg2 '2 2 xg2 zg2 

+ K 6 + D 0 - T_0 

y y y y d x 

Now, substitute (2-8), (2-48) and (2-50) into (2-54) to get the rotor 
moment equation in terms of shaft velocities 


M « A(to +0 - N0 ) + (C - B) (to +0 + N0 ) N + 

x x x y y y x 

A. (to +0 /2) + (C. - B_ ) (u + N0 /2) N + 

1 x x 1 1 y x 

A 0 (to - N0' +0 /2) + (CL - B ) (io + N0 /2 + 0„) H + 

2 x y x , 2 2 y x y 


K , * 0 /2 + D • 6 /2 + T_0 
xl x xl x D y 


M = B(to + 0 + ue ) - (c - A) (to +0 - N0 ) + 

y y y x 7 x x y 

* * n 

B 0 (co + N0 /2 + 0 ) - (CL - A.) • 

2 y x y 2 2 

(to - N9 + 0 /2) N + K 0 + D 0 - 

x y x' 7 yy yy Dy 


(2-56) 


Collecting terms in ( 2 - 56 ) gives 
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M x = ®x EA + (A i + A 2 )/2] + 6 x V 2 


e x Or[c - B + - B 1 + C 2 - B 2 )l/2] + K^/2} 

9 y N(C - A - B + C 2 - A 2 - Bg) + 0 y T D + 

to (A + A + A_) + to N[C - B + C, - B + CL - B ] 
x 1-2 y 1122 


M = 9 [B + Bj + 0 D +8 {N (C -A + C„-A_) + K}- 

yy 2 y y y 22 y 

0 x K(C - A - B + Cg/2 - Ag/2 - Bg/2) - 
6 x T D + W y (B + V " (C - A + Cg - Ag), 


7 (2-57) 


which are the coupled differential equations of motion in shaft fixed 
coordinates for the rotor of a two-gimbal gyro. 

Complex Transformation 

Proceeding as before, the symmetrical components are defined 


I = [A + B + (A^ + Ag)/2 + Bg]/2 
I- = [A - B + (A 1 + Ag)/2 - Bg]/2 


D a [D _/2 + B ]/2 
xl' y 

IT a [D _/2 - D ]/2 
xl y 


K - [K^/2 + Kj/2 

K- " [K _/2 - K 3/2 
xl y 

J 53 (2C - B - A 3Cg/2 + C /2 - 
J- = [ - b - A - (B 1 + Bg + C 1 


(2-58) 

\/2 - Bg/2 - Ag)/2 
+ Cg)/2 + Ag]/2 I 


L = [2(C - A - B) + 3(Cg - Ag - Bg)/2]/2 

L' = (Cg - Ag - Bg) A 
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P - (A + B + + Ag + Bg)/2 

P - = (A - B + A ± + A 2 - B 2 )/2 

R = (2C - A - B + 2Cg + ^ - B^ - B 2 - Ag)/2 
R* = (A - B + C ± - B 1 - B 2 + A 2 )/2 

These symmetrical components (2-58) are substituted into (2-57) to 
give the rotor moment equation in terms of symmetrical components, 
which is identical symbolically to (2-2l), only the components are 
defined differently in (2-58). As in the single gimbal case, (2-21 ) 
is combined into a single equation using complex notation, transformed 
to the case-fixed coordinate system, and Laplace transformed to give 
(2-28). The coefficients of variables, Z^s), Z 2 (s), etc,, are defined 
in (2-24) and substituted into (2-28) to give 


Z 1 (s)0 XY (s) + Z 2 (s)0 XY (s - J2H) = Mj^fs) - 

Z 3 (s)<j. XY (s) - Z 1 + (sH xy (s - J2R), 


(2-30) 


which is manipulated to give a conjugate equation at twice spin fre- 
quency. 


Z 1 (s - j2R)0 XY (s - J2N) + Z 2 (s - j2N)0 XY (s) = M^s - j2N) 

-Z 3 (s - j2N)? XY (s - J2N) - Z 4 (s - j2N)<j. xy (s> . (2-31) 

Simultaneous solution of (2-30)and (2-31) yield 


©■^(s) 

0^(5 - J2N) 
i 3 (s) 

Z u (s - J2H) 


-1 


\(s) 

Z S (B) 


Z 2 (s - J2H) 

Z x (s - j2N) 

Z 4 (s) 



Z 3 (s - j2N) 




V <s) 


M^ts - j2K) 


(2-32) 



kl 


Again the simplifying assumptions of 

1. Negligable drag and. damping, therefore 

T = D = D' = D = 0 

V K 

2. Symmetry, A = B, A^ = B^ , - B^, and 

K=K=K=k 
•xl x2 y 

allow us to evaluate (2-32) by simplifying (2-58) which becomes 


I = A + AjA + 3 A A 
I" = A 1 /U - A 2 /k 

K = 3 A k 

K" = - i/k k 

J = c - A + 3 C A - C 1 /U - Aj/4 - 3 A 2 /4 
* (C x - C 2 - A 1 * A 2 )A 

L = C- 2A + 3 C 2 A - 3 A 2 /2 
L“* = C 2 /4 - A 2 /2 

P « A + A 1 /2 + A 2 
P' = A^/2 

R=C-A+C 2 + C 1 /2 - A 1 /2 - A 2 
S' = C '/2 - A 1 /2 


S (2-59) 


Which may in turn be used in (2-29) to obtain 
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Z 1 (s) « s 2 (A + Aj/ 4 + 3A 2 /4) - j2tts (C/2 + 3C 2 /8 

+ a^/ 4) +(3/4)k- n 2 (a ± /2 - c 1 /4) 

Z 2 (s) = s 2 (A 1 - A 2 )/4 - jNs(A 1 /2 - C 2 /4) - k/4 
- n 2 (a 1 /2 - c 1 /4) 

Z 3 (s) - s(A + j^/2 + Ag) - JH(C + C x /2 + Cgj 
Z k (s) = sA 1 /2 - ~ C x /2) . 


> (2-60) 


Conjugation and substitution of s - j2N for s in (2-60) gives 


Z^s 


- j2N) = s (A, + k ± /b + 3Ag/4) - j2Ws(2A - C/2 


'I 


+ A =l /4 + 3A 2 /2 - 3Cg/8) +&/V)k - N [4a - 2C + A 1 /2 

(2-61) 


+ 3A 2 - 3C 2 /2 - C 1 /4] 


Z 2 (s - j2N) = s (A^ - A 2 )/4 - jNs(A 1 /2 - A £ + C g /4) - k/4 
- N 2 (A 1 /2 - Ag - C 1 /4 + C g /2) 


Z 3 (s - j2K) = s(A + A 1 /2 + Ag) - jH(2A - C + A^_ - Vj2 - C g ) 
Z 4 (s - J2H) = sA 1 /2 - JITC 1 /2 


Transfer Function 

The transfer function for the two gimhal gyro may be found by 

using 


e (s) « - jafjM^s) - Z 2 (s)^y.Cs - J2H) 

~[¥ 1 (s - J2H)Z (s) - Z 2 (s)Z }( (s - j2N)]<j^(s) 

-[^(s - j2K)Z^(s) - Z 2 (s)Z 3 (s - j2N)]^y(s - j2N)> 
/CV s) • Z 1 (s - J2H) - Z 2 (s) • Z 2 (s - j2N)3, (2-36) 



along with (2-60) and (2-6l). The approximation used is that rotor 
inertias are large compared with gimbal inertias. The constant term 
of the characteristic equation, Z^(s) • Z^(s - j2N) - Z^(s) • Zg(s - j2W) 
is 


[3/4k - - Oj/4)] ♦ [3/Uk - 

- 3C 2 /2 - c x A3 - t-k/4 - m 2 (a 1 /2 

- w 2 (a 1 /2 - Gj/4 - a 2 + c 2 A)] 


N 2 (4A - 2C + A 1 /2 + 3A 2 

- C/4)] . [-k/4 

1 (2-63) 


which is forced to zero. It is expedient to do this in the following 
manner : 

1. Make the constant term of Z (s) equal zero. This 
removes rotor inertias fronithe tuning condition, 
allowing the designer independent choice or rotor 
and suspension characteristics. 

2. Make the constant term of Z (s - j2N) equal zero, 
completing the tuning procedure. 

Thus 


and 


3 /4k - K 2 (A 1 /2 - CLj/4) = 0, 


k/4 - H 2 (A 2 - C 2 /2 - + 0^4) *» 0, 


> (2-64) 


which may he realized if 


A 2 " C 2 /2 " 2A 1 /3 + C l /3 “ 0 
one solution of which is 


(2-65) 


A 2 = 2/3A 1 , C 2 = 2/3C 1 


( 2 - 66 ) 



kb 


Equation (2-36) may now be written 


6 xy (i!) ' V U> 


z a (6) 




z x (s - 4211)^(3) 


M (s - J2W) 
XY 


(2-6T) 


Zo(s) : 


*^(8) 

L Z ( s ) 


zTGT'** 


Z 4 (s) 


Z £ (s) Z 3 (s - j2N) ^ ^ 

~ Z~( s) Z 1 (s - j2IT) *XY^ S " 


as 'the timed transfer function. .Under this condition 


Z 1 (s) « b A - j2NsA 

Z x (s) • Z (s - 32N) = A 2 s 2 (s - j2N) 2 

Z 3 (s) = sA‘ - j2NA 

Z 3 (s)/Z 1 (s) = 1/s 

Z 3 (s - j2ST)/Z 1 (s - j2N) = l/(s - J2H) 


and the transfer function becomes 


V s) = 




s) 


sA(s - J2N) 


A 1 /12 M (s - J2N) 
A 2 s(s'->j2ff) 


'<>- 68 ) 


^(s) 5 \ ^(s - J2K) 


M. 

s 


12 


A (s - j2Ef) 


Separating (2-68) into its real and imaginary component 


0 x (s> = !V!2 


2R 


• My(s) 




As(s P + 4W 2 ) As(s 2 + kM 2 ) 


A^/12 sM^s - j2If) + 2ff • M y (s - j2M) 


s(s 2 * UH 2 ) 


5 A 1 S< ^ S ” ^ 21 ^ + 2N * ” ^ 21 ^ 

12 . . 2 * „2 v 

A (s + 4N ) 


(2-69) 



1*5 


2N * Mj^(s) sMjj(s) <j)y( s ) 

-J- .... . - .... ^ — ■ ■■ 

As(s 2 + UN 2 ) As ( s 2 + 4N 2 ) s 

2MM x (s - j2N) - sMjj ( s - J2N) 
s(s^ + 4-N 2 ) 

2N ♦ t x (s - J2N) - s? y (s - J2H) 

(s 2 + 4 n 2 ) _ 

The transfer function for direct rate inputs , and moment inputs , hoth 
direct and 2N is similar to that of the single gimbal gyro (2-44 ) , . 
differing only in the scale factor of the 25? moment input. A constant 
angular displacement at 2N (impulse of rate) is not rectified, which 
improves the quality of the gyro as a rate integrating sensor. . - 

V. COMPARISON OF ONE AND TOO GIMBAL GYRO DYNAMICS 

The transfer function for direct rate and moment inputs, which 
* 

is required by the system designer to use the gyro as a strapdown 
instrument, is identical for both the single and two gimbal gyros. 

The principle differences are (l) elimination of the undesirable 2N 
rectification effect, and (2) a tuning condition, (2-64), which is 
independent of rotor inertias, in the two gimbal gyro. Both of these 
characteristics are important. The 2N rectification in a single gimbal 
gyro destroys its usefulness as a high quality attitude sensor. A 
tuning condition which contains rotor inertias (2-Uo) means that the 
dynamic antispring is not a function of gimbal inertias and speed 
alone, consequently, the rotor cannot be considered free, the re- 
quirement of a high quality rate integrating gyro. This lack of 
freedom is displayed in the 2N rectification of the single gimbal gyro. 
The conclusion which may be reached here is that a tuned gyro must 
have more than one gimbal for successful operation as a sensor. 


V 12 




CHAPTER III 


REBALANCE LOOP ANALYSIS FOR A DRY, TUNED, TDF GYRO 

For strapdown applications, the plant dynamics derived for the 
•gyro of the preceeding chapter must he included in a closed loop. This 
loop causes the gyro to be a null instrument, any sensed deviation of 
the rotor position with respect to the ca'se from null causes an applied 
restoring or rebalance torque to force the rotor bach to null. As in 
the case of the SDF gyro, the torquer current can be measured to give 
an indication of sensor motion. The plant dynamics, including torquer 
coil dynamics and pickoff gain, is shown in Figure 3-1. Both analog 
and the digital rebalance loops will be investigated. 

I. ANALOG LOOPS 

For ease of design and hardware implementation, rebalance loops 
using continuous signals are preferred. Two different analog rebalance 
loop designs are presented in this section. 

Teledyne Loop 

The Teledyne design is based on meeting the following criteria: 

1. The steady-state output angle (rotor hangoff) errors ' 
are zero for constant angular velocity inputs . 

2. The maximum absolute transient error is less than 5 

. milliradians to keep the rotor from hitting its stops. 

3. The feedback loop gain must be attenuated at the spin 
frequency, N, to avoid rectification errors at this 
frequency. 

This design was implemented using the plant transfer function in complex 
form, 

0 , , "ey 131 W s) 

0 XY^ As ( s - J2H) 
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(2-38) 
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A root locus design was used resulting in the rebalance of Figure 3-2, 
a single loop with the following variables in complex form: 

Rate input <j>-jry(s) 

Moment input M X y(s) 

Output (rotor hangoff) angle 0yy 


This loop contains physically unrealizable poles and zeros, that is, 
complex roots which are not accompanied by a conjugate. Using 2-22, 
which converts the variables in complex form into their real variable 
components , the feedback component of the loop of Figure 3-2 is trans- 
formed to Figure 3-3 in which both direct and cross-coupled terms 
appear. The general effect of complex roots without conjugates in 
the complex variable domain is to produce both direct and cross- 
coupled terms in the real -variable domain. 

Honinter acting Loop 

Another analog rebalance loop uses the well known principle of 
Q 

noninteraction, in which the matrix of the open loop transfer function 
is diagonalized. This diagonalization causes the closed loop matrix 
to also be diagonal, which decouples the response of all outputs but 
one to a given input . Thus , each, input is paired with an output , and 
these input-output pairs are noninteracting with each other. This 
noninteraction essentially reduces a multi-variable system with n 
inputs and n outputs to n separate single-input-single-output systems 
which may be compensated individually using classical techniques . 

The simplest method of diagonalizing a matrix, M, is to multiply 
it by its inverse, M -1 , since MM -1 = I, provided M is non-singular. 

The plant gain matrix (which excludes the torquers and pickoffs) is 


G = 


1/A 


-2N/s 


s 2 + 4 n 2 


2N/e 


(3-1) 




Figure 3-2. Teledyne Rebalance Loop in Complex 

Coordinate Form 4 =- 








Rebalance Loop in Realizable Form 
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which, inverted, "becomes 


(3-2) 

Because of linearity of the system, this inverted plant gain matrix may 

be placed anywhere within the feedback loop to obtain diagonalization. 

2 

In addition, each element of this matrix may be divided by s, or s , 
to give a type I, or II, closed loop response, respectively. 

In order to compare the performance of the noninteracting loop 
with that of the Teledyne loop a type II system was designed. The 

C 

Bode plot of the open loop is shown in Figure 3-4. The gain of ICr 

gives a steady-state error in the rotor hangoff angle of 0.5 milli- 

' 2 
radian for an acceleration input of 50 rad/sec , which is comparable 

to the Teledyne response. A lead network is used to add a zero at 

a) = 200 so the slope of the open loop transfer function is -20 db/ 

octave at the zero db crossing. This gives a theoretical 10 to 90 

percent rise time of 7-3 msec and a 1.2# overshoot. The additional 

poles in the transfer function are at to = 10000 rad for the torquer 

coil and w = 20000 rad for the lead network pole. This design permits 

the addition of a notch filter with imaginary zeros a w = +_ j N, and 

two real poles at oj = 3N, 3SF = 628 rad,- if one is considered necessary. 

The single-input-single-output transfer function is 


= As 9 


2N 


-2N 


6 x 10 5 (s/200 + l) /. 

• “ O * ' 

<f> s (s/20000 + l) (s/10000 + l) 

X 

Figure 3-5 shows a block diagram of the noninteracting loop used to 
achieve this transfer function. 






Figure 3-5. " Uoninteracting Analog Rebalance Loop 
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Comparison of the Teledyne and Noninteracting Loops and the Results 
of Simulation 

Both the Teledyne and noninteracting designs are satisfactory to 
rebalance the dry tuned TDF gyro. The Teledyne loop, -which in realizable 
form is a two-input-two-output system, is the more difficult of the design 
methods to implement. The implementation of this design involves using 
system variables in the complex coordinate form (Figure 3 - 2 , page ^9 ) 
which results in a root locus design in which real-axis symmetry does 
not exist. The notch filter is an essential part of this design; its re- 
moval causes the loop to become unstable. 

The noninteracting design, which reduces the design problem to 
that of a single-input-single-output system, has one flaw. The non- 
interaction depends on exact knowledge of decoupling and plant param- 
eters, in real world problems this never the case. The consequences 
of this inexact noninteraction are imaginary axis closed loop 'poles 
in the vicinity of 2N, and small off-diagonal elements in the transfer 
matrix, which also contain imaginary closed loop poles near 2N. These 
poles limit the closed loop bandwidth of the system. 

Results of simulation shown in Table 3-1 of the Teledyne loop and 
the noninteracting loop with a notch filter included shows the Teledyne 
loop having slightly superior performance. 

Table 3-1: Rotor Hangoff Angle for Dry £uned TDG Gyro with Acceleration 

Input (x-axis) of 50 rad/sec 


Rebalance Loop 
Loop 

Rise time 
10 %- 90 % 
msec 

5% Settling 
Time 
msec 

% Overshoot 

S . S . Error 
Millirad 

Honinteracting 

6.0 

21.0 

19.2 

0.50 

Teledyne 

6. 0 

15.0 

9.2 

0.53 
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II . DIGITAL REBALANCE 

Pulse width modulators (PWM's) may he inserted ahead of the 
torquers so that the loop may operate in .the PTSA mode. From the 
equal area principle"^, a linear system gives- the same response at 
the time of sampling, regardless of the waveshape of the input, pro- 
vided the area under the input curve is equal. Thus the current to 
the torquer may be pulse width modulated without a deterioration of 
sensor performance. 

One type of PWM employed is binary width modulation (BWM) , in 
which the torquer current assumes a constant magnitude of positive or 
negative polarity, such that the net current time produce (area) over 
a sampling period is proportional to the input of the BWM at the 
sampling time. Because constant current is delivered to the torquer 
coil, its power is constant regardless of the net current; consequently, 
unstationary thermal gradients in the torquer coils, which can cause' 
scale factor change, are avoided. Data resolution is obtained by 
varying the positive and negative widths in small discrete steps , as 
.showru in Figure 3-6 . 



H h Resolution of Current 

— n ' r- 


Sampling 

Period 



Figure 3-6. BWM Current 
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Both the Teledyne and noninteracting designs were simulated with 

BWM's in the loop ahead of the torquer coils (Figures 3-7 and 3-8). 

The results of these simulations are shown in Figure 3-9 through 3-12 

for sampling periods of 0.5 and 0.1 msec. A constant angular accele- 

2 

ration input about the x-axis of 50 rad/sec is used in each case. The 
x- and y-rotor hangoff angles and inputs to the x- and y-BWM's are 
plotted. 

For a sampling" period of 0. ‘5 'msec, the Teledyne loop produced 
objectionably large 500 Hz ripple at the inputs to the BWM's , as well 
as substantially large transients in the rotor hangoff angles. For 
the noninteracting loop the inputs to be- BWM's are' free of ripple, 
but have constant offsets from the analog case. The rotor hangoff 
angles have approximately 4 mrad peak-to-peak ripple at 200 Hz. 

Therefore, binary width modulation with a sampling period of 0.5 msec 
is considered unsatisfactory for both the Teledyne and noninteracting 
designs . 

For a sampling period of 0.1 msec- the Teledyne loop has a small 
amplitude disturbance of undiscernable low frequency in both the rotor 
hangoff angles and the inputs to the BWM's. The variables approach the 
values of those of the case of analog rebalance, and the small amount 
of deviation is not considered sufficient tp degrade, sensor- performance . 

The noninteracting loop has a 200 Hz frequency component in the rotor 
hangoff angles of approximately 0.4 mrad peak-to-peak, which is acceptable. 
The BWM inputs contain no noticeable deviation from the analog case. , For 
both the Teledyne and noninteracting designs, binary width modulation 
with a sampling period of 0.1 msec is acceptable. 

Existence of Limit Cycles 

In designing a PTSA loop, it is desirable to avoid limit cycles, 
since they can cause signal excursions beyond acceptable limits of 
operation of the system and can also generate extraneous information 
regarding sensor motion. This section is devoted to the determination 
of the existence of limit cycles using a describing function (DF) for 
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the non-linear element in the loop, in this case the PWM. ■ DF theory is 
well documented in the literature 11 for single non-linearities within 
a system, and involves an approximation which linearizes the non-linearity 
Ordinarily, a single frequency input to the non-linearity produces an 
output including this frequency and its harmonics,. Most servo loops are 
'considered to he good low pass filters, and the higher harmonics are- 
assumed to he attenuate hy the loop. Consequently, only the • fundamental' 
frequency of the output of the non-linearity is considered, and its DF 
may he defined as 


N(A,fj)) 


Fundamental Harmonic of Output 
<■ Input 


(3-4) 


The system of Figure 3—13 is separated into. its linear L(to), and non- 
linear, N(A,<j> ) , components where ‘A is input amplitude t-0 the non—’ ■ . 
linearity," to 'the -frequency, 'and <j> the phase •angle - of N, and if the input 
to the system is zero, 

: 1 + L('w) N(A>) = 0. (3-5) 


External 

Input 




Linear Components 
L(to) 




I'Ton-Linear Components , 
N(A,4>) 


Figure 3-13. Equivalent Loop With Non-Linearity 
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Equation (3-5) is generally solved graphically by plotting the linear 
portion L(to) on an amplitude-phase plot of -l/N(A,<f»). Permissible .limit 
cycles occur at the intersection of the L(co) and -l/N(A,<f>) curves. 

In a sampled system, further constraints are placed on the per- 
missible limit cycles. The period of the limit cycle must be some in- ‘ 
t'egral multiple of the sampling period, and in the case of a loop which , 
contains at least one integration, this multiple must be even. 

For a sampled system using a BWM as its non-linearity. Figure 3-1 4 
shows the lower bounds of limit cycles wi-th periods of 2 through 8 times 
the sampling period. If the operating point of the linear -portion-, L(w), 
lies above one of these boundaries, thep that limit cycle can exist; 

For systems with multiple non-linearities, 'DF theory generally 
cannot be used, unless some unusual property is present. Fortunately , 
the noninteracting design effectively decouples one BWM from the other ,, 
so that DF theory may be used on the single-input-single-output equiva- 
lent loop. Figure 3-15 shows a plot of (3-3) on Figure 3-l4. Because 
(3-3) represents a type II system, its phase approaches -l80° for low 
frequencies, making low frequency limit cycles unavoidable. For this 
reason a type I' system is desirable, and, modification of the feedback 
loop to make the system type I includes a lag network, w-ith a pole at . 
s = -14 and zero at,s = —200, so that the closed .loop bandwidth can 
remain unchanged. The transfer function for this type I system is 


x K } 8000(s/200 + l) 

* / . s(s/l4 + l)(s/l0000 + l) 

Ms) 


(3-6) 


and is also plotted in Figure 3-15* At low frequencies, the phase of 
this type I system approaches -90°, and low frequency limit cycles can 
be avoided. For a digital rebalance loop employing BWM,, the linear 
portion of the loop must be type I. 
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Comparison of Types I and II Loops 

In practice, the requirement of zero steady-state rotor hangoff 
error to constant angular velocity input is unnecessary, a small con- 
stant error can he tolerated. Although angular acceleration inputs are 
experienced in practice, their time integral is sufficiently small that 
the rotor -will not hit its stops ifr a properly -designed loop. 

Regardless of whether the response of the rotor hangoff angle 
to case motion is type I or II, the response at the input to the 
torquers is the same for a given motion input. The principle of con- 
servation of momentum causes the gyro to process with an angular 
velocity, co, when torque, T, is applied, according to the relation 


T = m x H (3-7) 

where E_, the rotor momentum, is the product of the spin velocity and 
spin-axis moment of inertia. For a constant velocity input to the 
x-axis, a constant restoring torque from the y-axis torquer is neces- 
sary to rebalance the gyro, 'independent of whether steady-state rotor 
hangoff angles are allowed for this input. 



CHAPTER IV 


PROPOSED PTSA LOOP FOR THE DRY TOTED TDF GYRO AND 
CONCLUDING REMARKS 

‘ T. DIGITAL REBALANCE IMPLEMENTED FOR THE TELEDYNE GYRO 

•The digital rebalance loop of Figure 4-1 gives the single-input- 
single-output transfer function. 


V ' _ ' 8000 (s /200 + l) 

: / v s(s/l4 + 1) (s/10000 + l)(s/20000 + l) 

( M S ) 


(4-i) 


which is similar to (3-6) except for the pole at s = -20000. This 
response is type I, but the pole at s = -l4 gives the same mid- 
frequency response as the type II system described by 3-3. 

The Teledyne gyro has the physical characteristics shown in 
Table 4-1. 


Symbol 

K 

P 

A 

’ C 

N 

H 


Description ' 

Torquer Scale Factor 
Piekoff Scale Factor 
Cross-axis Moment of Inertia 
Spin-axis Moment of. Inertia 
Spin Frequency 
Angular Momentum 
Torquer Coil Resistance 


Value 

20.1° /amp. sec 
2.26 volts/ 0 
8 00 gm cm^ 

2 

1500 gm cm 

100 Hz (638 rad/sec) 

. J> 2, 

10 gm cm /sec 
12 ohm 


Table 4-1: Teledyne Gyro Characteristics. 
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These characteristics constrain the digital rebalance loop, -which is 

12 

the U.T. modification of the Hamilton Standard design. This loop, a 
BWM design originally developed for floated, gas gpin "bearing, SDF gyros, 
has been adapted to accelerometers as well. 

From the physical characteristics of the gyro, the non-interacting 
design transfer' function of (1— l) and the' design procedures outlined in 
Reference 12 for the U.T. loop, the following considerations must be 
observed: 

1. Two channels of the equivalent rebalance loop for an 
SDF gyro are needed, as well as additional components 
for decoupling. 

& 2 

2. The rather large rotor angular momentum of 10 gm cm /sec 
requires a large rebalance torque. A rate input of 6o°/sec 
requires 3 amperes of torquing current and over 100 watts 
of power to the torquer coil. 

3- The ramp slope constraint simply assumes that the error 
signal at the input to the BWM does not change faster 
than the ramp with which it is compared. Thus the ramp 


R > K K, K I, " 
pc t p 


(13 of Ref. 12) 


where K does not include the gain A in the decoupler.. 

4. The interrogation or sampling period, t., of the BWM 
puts a more stringent requirement on the closed loop 
bandwidth via the ramp slope constraint than through 
the Nyquist sampling frequency. This relation is 

f . 

f - = 7^ (16 of Ref. 12") 

Cl 7T 

max 


5 . 


where f. is l/t.. 
i i 

The gain, K , is found from (8) of Ref. 12 at u = 1. 

po 


G(io = 1) 


2 K, K K I 
t pc p 


Rt. 


l 

where I is the torquer current in amperes, 

8000. K does not include the gain A, 

PC 


(U-2) 

and G(m = l) is 



71 


Figure 4-2 shows a block diagram which uses the U.T. rebalance 
loop with a few additional components. The gain of 24-000 is the produc 
of A = 800 and K = 30 from (4-2). The interrogation period, t^, is 
0.1 msec to avoid limit cycles and provide sufficient closed loop band- 


width . 


The data resolution of the loop can be stated as 


Ae_, = 7200 
d 


Maximum Rebalance Rate 
Data Rate 


(arc sec) 


(4-3) 


and is 0.7 arc sec for a 60°/sec maximum rebalance rate and 6l4.4 kHz 
data rate. This resolution may be improved by dual mode operation of 
the torquer current, or by increasing the data rate. 


II. CONCLUSIONS 

This report includes the investigation of the feasibility of 
using a JPTSA rebalance loop for a dry, tuned TDF gyro in the strapdown 
mode. This strap do™ TDF gyro would be a welcome addition to the famil 
of sensors available for spacecraft navigation. 

Also included is a derivation for transfer function of this gyro 
including the conditions of tuning for both the single and two gimbal 
cases. The single gimbal gyro is not a true "tuned" gyro, since its' 
tuning condition contains rotor inertia terms,' except in the special 
(and physically unrealizable) case where the rotor is infinitesimally 
thin, causing the C-2A term involving rotor inertias to vanish' from 
the tuning equation. In the strict sense, tuning is a function of 
speed, gimbal inertias, and the torsional spring constants of the 
restraints, and this can occur for a physically realizable rotor only 
when the gyro has more than one gimbal. 

Two analog rebalance loops for this gyro were simulated. These 
loops were converted to the digital mode by the addition of binary 

t 

width modulators to digitize the torquer currents. Although these 
loops performed well at the higher sampling rate of 10 kHz, the 
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recommended final design for a digital rebalance loop for this gyro 
should be Type I, to assure- operation without low frequency limit 
cycles. 

A conceptual type I design using a modified version of the U.T, 
rebalance electronics is given in the preeeeding section of this chapter. 
While this design is feasible the gyro used in -the loop has poor charac- 
teristics for a high rate environment, requiring excessive power to the , 
torquer coils. However, the concept of digital rebalance for a dry tuned 
TDF gyro is sound, and awaits a gyro with suitable characteristics. 
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A USEFUL COMPLEX CONJUGATE PROPERTY 

If B{s) - A{ s - jy), then B(s) = A(s + jy). 

Proof : 

Let j he the imaginary rotation for the spatial domain, and i he the 
imaginary rotation for the time domain. The Laplace transform variable 
s = y + ig. 


a(t ) = x(t) + jy(t), A(s ) = X(s) + jY(s) 
therefore 

a(t ) = x(t) - jy(t), and A(s) = X(s) - jY(s). 


Let 


h(t) = a(t)e^ Yt = (x( t ) + jy(t))*e^ Yt 

then 


h(t ) = a(t)e“^ Yt = (x(t) - jy(t))*e“^ Yt 


and 


B(s) = X(s - jy) + JY(s - jy) = A(s - jy), 
B(s) « X(s + jy) - JY(s + jy) = A(s + jy). 


Q.E.D. 



